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ABSTRACT 

We estimate the accretion rates of 235 Classical T Tauri star (CTTS) candidates in 
the Lagoon Nebula using ugriRa photometry from the VPHAS+ survey. Our sample 
consists of stars displaying Ha-excess, the intensity of which is used to derive accretion 
rates. For a subset of 87 stars, the intensity of the u-band excess is also used to estimate 
accretion rates. We find the mean variation in accretion rates measured using Ha and 
w-band intensities to be ~ 0.17 dex, agreeing with previous estimates (0.04-0.4 dex) but 
for a much larger sample. The spatial distribution of CTTS align with the location 
of protostars and molecular gas suggesting that they retain an imprint of the natal 
gas fragmentation process. Strong accretors are concentrated spatially, while weak 
accretors are more distributed. Our results do not support the sequential star forming 
processes suggested in the literature. 

Key words: accretion, accretion discs, stars: pre-main sequence, stars: variables: 
T Tauri, Herbig Ae/Be, open clusters and associations: individual: NGC6530, Lagoon 
Nebula, M8 


1 INTRODUCTION 

Classical T-Tauri stars (CTTS) are low-mass pre-main se¬ 
quence (PMS) stars thought to accrete mass from a circum- 
stellar disc via a stellar magnetosphere. Mass accretion rates 
(Macc) are thought to decrease with time as the circumstel- 
lar disc material is depleted. Simple viscous disc evolution 


predicts the observed trend rather well I Sicilia-Aguilar et al. 
|2006| |Barentsen et al7||2011| |Manara et al.||2012| |. Obser¬ 
vations also reveal that M acc scales steeply to the square 


of the stellar mass M* l 

Muzerolle et al. 

2003j Natta, Testi 

& Randich 2006. Sicilia-Aguilar et al. 2006. Rigliaco et al. 

2011 

Manara et al. 2012 

Kalari & Vink 2015 1 in a variety 


of environments, however no such dependence is actually 
predicted by current theory (H artmann et ah] [2006). Ob¬ 
servational M acc estimates can provide an empirical base 
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against which disc evolution models can be tested, and star 
formation histories of individual star-forming regions can be 
de-constructed ( |Barentsen et al7||2011[ |Venuti et al. ||2014]) . 
Selection effects (Clarke & Pringle 20061 and intrinsic vari¬ 
ability ( Scholz et al.||2005 Costigan et al.|[2012 2014 1 may 
however affect the observed relations, and for these reasons 
large, uniform statistical samples are essential to constrain 
theory. 

Empirical M acc values have been derived from 


intermediate-to-high resolution spectra (Muzerolle et al. 


2003| Natta et al.|2006; Sicilia-Aguilar et al.|2006i), (7-band 

1 Rigliaco et al. 

2011 

1 or Ha photometry (Barentsen et al. 

2011 

De Marchi et al. 2010; Spezzi et al. 2012; Manara et al. 


2012). Statistical studies of M acc have generally been re¬ 


stricted to relatively nearby star-forming regions, such that 
the effects of either massive stars or the variation of the 
star-forming environment on PMS evolution remain poorly 
constrained. In particular, whether photoevaporation from 





















































2 V. M. Kalari et al. 


OB stars in the local environment affects M acc remains an 
open question ( Williams &; Cieza|2011 |. 

The advent of deep, high-resolution wide field narrow 
and multi-band surveys of the Galactic plane affords new op¬ 
portunities to study star forming regions in greater details 
and to uncover significant new samples of CTTS in a vari¬ 
ety of star-forming regions to address the problems raised in 
the preceding paragraphs. The VST Photometric Ha sur¬ 
vey of the southern Galactic plane and bulge (VPHAS+) 
is an ugr iB a ima ging survey reaching down to r ~ 21 mag 
( Drew et al.|2014 |. M a cc rates can potentially be estimated 
photometrically from either the u-band or Ha excess with 
known detection limits. It thus becomes possible to con¬ 
struct a sample of M acc rates that is relatively complete for 
PMS masses down to 0.5 Mg at 2 kpc, assuming moderate 
extinction (Ay ^ 2). This enables a uniform M acc data-set 
in a variety of notable star-forming regions, including the 
Carina Nebula, NGC6611, and NGC3603. Such work has 


already begun by Vink et al. (2008) in the Cygnus OB 2 


massive star forming region, Barentsen et al. (20111 in the 


IC 1396 star-forming region, and Barentsen et al. (2013) in 


the open cluster NGC 2264 using data from the counterpart 
northern Galactic plane survey IPHAS (D rew et al.|2005 |. 

In this paper, we aim to measure M acc for a signifi¬ 
cant fraction of CTTS candidates in the Lagoon Nebula 
(M8) identified using VPHAS+ survey photometry. M8 is 
a southern Hu region located in the Sagittarius-Carina arm. 
It comprises the embedded open cluster NGC 6530 located 
to the east of the central dust lane, and the optically ob¬ 
scured Hourglass nebula cluster to its west (see Fig. 1). It 


contains at least sixty known PMS members (Sung, Chun 
fe Bessell| 2000| IPrisinza.no et aH2005| |Arias, Barba fe Mor- 


rell||2007 |Tothill et al.|2008l, but with very few M acc esti¬ 


mates in literature (Gallardo, del Valle & Ruiz 2012J). The 
median age of PMS members in the region has been es¬ 
timated ~ 1 Myr using colour-magnitude diagrams (CMD; 


Sung et al. 

2000; Prisinzano et al. 2005| Mayne et al. 2007; 

Arias et al. 

2007), concurrent with both the dynamical age 

(Chen, de Grijs & Zliao|2007 van Altena & Jones|19721 and 


the main-sequence position of the confirmed O-type mem¬ 
bers ( |Toth ill et al.1 |2008| ). It is comparatively younger than 
most regions for which a statistically significa nt number of 
M acc estimates are known (e.g. see Fig. 8 in Spezzi et al. 


20121 to enable a meaningful statistical analysis. However, 


it has been suggested that a previous burst of star forma¬ 
tion occurred around ~ 10 Myr ago by Lada et al. (19761 


and van den Ancker et al. (19971. Lightfoot et al. (19841; 


Damiani et al. (20041; Prisinzano et al. (20051; Arias et al. 


(20071; Prisinzano et al. (20121 have reported sequential age 


gradients of the PMS members in region. 

The mean reddening E(B — V) towards NGC 6530 mea- 


fitting (van den Ancker et al. 1997 

Arias et al. 2006) and 

isochrone fitting (Prisinzano et al. 

2005 Mayne & Nay- 


these studies is ~ 0.5 mag for a few stars, and the fore¬ 
ground extinction ~ 0.2 mag. Distance ( d ) determinations 
to NGC 6530 using isochrone fitting provide d = 1250±50 pc 
(Prisinzano et ak1|2005 Arias et al.||2007| Mayne & Naylor 


2008). Previous measurements indicate d as high as 1800 pc 
( |van den Ancker et al.||l997| |Sung et al.||2000| . This dis¬ 
crepancy is thought to arise because larger determinations 


rely on fitting only O and early B type stars, rather than 
the complete blue edge of the distribution of stars in the 


upper main sequence on the CMD. Prisinzano et al. (2005) 


argue that any background stars lying behind the nebula 
must be significantly reddened, moving them away from the 
main sequence. Based on the above argument, |Prisinzano| 


et al. (2005) suggest that the blue edge of the distribution 


of stars in the CMD more accurately defines the zero-age 
main sequence at the distance to the Lagoon Nebula. 


Most authors (Neckel & Chini 1981 

Sung et al. 2000 

van den Ancker et al. 1997 Arias et al. 

2002 2006) have 


found the reddening to follow a normal reddening law, 
Rv «3.1, but found that H36, an 07 V star located in the 


Hourglass nebula with Rv ranging from 4 - 8 (Johnson 1967 
|Sung et al.|2000| | . This anomaly is explained by the presence 
of larger silicate grains surrounding H36, suggesting selec¬ 
tive evaporation in its immediate environment leading to an 
anomalous Rv in the close vicinity ( |Hecht elTa l. |1982| . Sung 
et al. (2000) used the E(B — V)/E(V — I) colour-excess ra¬ 
tio and found a Rv range between 3.0-3.6 for OB stars in 
the Lagoon Nebula with a mean of 3.1, but higher values 
for PMS stars which were effected by blue excess emission, 
and H36. We note that Prisinzano et al. (2012) used the 
longer baseline afforded by the E(B — V)/E(V — K) ratio 
to find Rv ~ 5, inconsistent with standard reddening. How¬ 
ever, many of the stars used to derive this reddening law 
were PMS stars that are most likely affected by excess blue 
and infrared emission. Prisinzano et al. (2012) note that 
because of this, they were unable to constrain Rv with pre¬ 
cision. We note that an anomalous reddening law of Rv = 5 
leads to a difference A r of s=s 0.5, which is approximately the 
same as the spread in E(B — V) values. 


At the mean reddening and distance to the Lagoon 
Nebula 0.3 Mq- 2 Mg PMS stars span a magnitude range 
13<r<21, making VPHAS+ photometry suited to esti¬ 
mate their M acc . The four confirmed O-type members plus 
some sixty B-type stars (Tothill et al. 20081 represent a step 
toward increasing the massive star content compared to sim¬ 


ilar works using IPHAS in the Northern Galactic plane (Bar- 
entsen et al.|20lT 20131. 


In Section 2, VPHAS+ observations and source selec¬ 
tion criteria are detailed. In Section 3, we identify CTTS 
candidates. We measure the accretion rates and stellar prop¬ 
erties of the CTTS candidates in Sect. 4. The distribution of 
M acc with respect to stellar mass is discussed in Section 
5. The spatial distribution of CTTS, and their properties is 
presented in Section 6. Finally, our conclusions are presented 
in Section 7. 


2 DATA 

VPHAS+ imaging is obtained by the OmegaCAM CCD im¬ 
ager mounted on the 2.6 metre VLT Survey Telescope (VST) 
on Cerro Paranal, Chile. The imager captures a 1° square 
field of view at a resolution of 0.21" pixel -1 . Each point¬ 
ing is supplemented by at least one offset exposure to min¬ 
imise CCD gaps. Imaging is carried out through broadband 
ugri filters and a purpose-built Ha filter. The central wave¬ 
length and bandpass of the Ha filter are 6588 and 107 A re¬ 
spectively. Exposure times are 150, 30, 25, 25, and 120s 
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Right Ascension (J2000) 


Figure 1 . VPHAS+ riHa three-colour image of the Lagoon Nebula. North is up and east is to the left. The Hu region Lagoon Nebula 
(or M8) is marked by the box having dimensions 50' X 30'. 


respectively. VPHAS+ reaches a 5cr depth at Ha = 20.5- 
21.0 mag and g = 22.2- 22.7 mag. Practical constraints have 
meant that the blue ( ug ) and red (riHa) observations are 
carried out separately. An additional r observation is car¬ 
ried out with every blue observation to serve as a linking 
reference. 

Observations used in this paper are detailed in Ta¬ 
ble [I] Data from the survey are pipeline processed at 
the Cambridge Astronomical Survey Unit (CASU), where 
nightly photometric calibrations are performed. Astrometry 
of sources is refined by matching with the 2MASS catalogue 
( Cutri et al.|2003 |. 

Overlaps between the observations were used to bring 
the red and blue data onto a common internal scale. 
The calibrations of band-merged catalogues were checked 
and refined using external gri data from the APASS 
all-sky catalogue (http://www.aavso.org/apass). We cross- 
matched bright sources (13 < mag < 15) morphologically 
classified as stellar in the VPHAS+ bands, and having 
random photometric error in each catalogue less than 
0.1 mag. Around 200 sources were selected for compari¬ 
son in each band. The median magnitude offsets found 
were Ar(red) =—0.027, Ar(blue) = —0.029, Ai = — 0.0106, 
Ag = 0.044. Aw = —0.23 mag was determined following the 
procedure described in Drew et al. (2014). 


2.1 Source selection 

The area of study is centred on RA 18 ,l 04 m 12 s , Dec 
-24°18'36" (J2000) and covers 50' x 30' (Fig. 1). This re¬ 
gion encloses the Lagoon Nebula (Tothill et al. 2008). We 
apply the following criteria to select sources: 


Table 1 . Log of observations 


Right Ascension 

Declination 

Filters 

Obs. Date 

Seeing 

18 h 02" 

a 06.32 s 

—24° 11/36.8" 

riHa 

2012 June 10 

0.6" 

18 fc 01" 

a 2SA7 s 

-24°00'36.0" 

riHa 

2012 June 10 

0.6" 

i8 h 06" 

a 18.95 s 

-24°11'36.8" 

riHa 

2012 June 27 

0.7" 

18 fe 05" 

a 36.10 s 

-24°00'36.0" 

riHa 

2012 June 27 

0.7" 

18 h 02" 

a 06.32 s 

—24°11'36.8" 

ugr 

2012 Aug 14 

00 

o 

18 fc 01" 

ri 2SA7 s 

-24°00'36.0" 

ugr 

2012 Aug 14 

00 

o 

18 fe 06" 

* 18.95 s 

—24°11'36.8" 

ugr 

2012 Aug 15 

0.9" 

is^os" 

* 36.10 s 

-24°00'36.0" 

ugr 

2012 Aug 15 

0.9" 


(i) r > 13 mag in both the red and blue filter sets to avoid 
saturated sources; 

(ii) random photometric error must be less than 0.1 mag 
in g,r (red and blue),i and Ha to keep photometric and prop¬ 
agated uncertainties small; 

(iii) morphologically classified as stellar in g and i and 


either stellar or star-like in r (sec Table 7 in Gonzalez-Solarcs 
et al.|2008 1 to avoid extended objects. 


(iv) morphologically classified as either stellar, star-like 
or extended in Ha and M-band photometry. The reasoning 
behind these criterion are given in the following paragraph. 

(v) average pixel confidence value greater than 90 and 95 
percent in broadband ugri and Ha respectively (see Drew 
et al. 2014). This is especially important in Ha to ensure 
the photometry remains unaffected by extra vignetting by 
the T bars separating the Ha filter segments. 


To filter for background variations, the CASU nebu- 
liser routine is performed on the VPHAS+ images as part 
of the data processing pipeline. The background is modelled 
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Figure 2. The g versus (g — i) CMD. The diagram includes 
all point sources meeting our selection criteria in g, i located in 
the Lagoon Nebula (see Fig.l). A reddening vector for Ay = 1 is 
shown. 



using a sliding median filter over a scale of 15". Although 
rapidly spatially varying nebulosity on shorter spatial scales 
can lead to imprecise sky subtractions. This is especially per¬ 
tinent when considering the variation of nebulosity on small 
scales as is the case in the regions surrounding the NGC 6530 
cluster and the Hourglass nebula (see Fig. 1). Imperfect sky 
subtractions or imperfect flat-fielding in such cases can af¬ 
fect the point spread function of sources especially in Ha? 
and u-band photometry respectively. In such cases, the fi¬ 
nal morphological classification of Ha or u-band sources as¬ 
signed by the pipeline may be +1 (extended), reflecting the 
resulting mild degradation of the point spread function. So 
as not to reject too many detections of sources that are 
otherwise clearly point-like, we relax the morphological con¬ 
straints and accept ‘extended’ Ha and u measurements. In¬ 
stead, we inspect the images of each selected Ha excess emis¬ 
sion line source by eye. An unsharp mask filter was applied 
to the reduced and calibrated single-band mosaics to sup¬ 
press large scale low frequency features such as diffuse back¬ 
ground, and highlight the high frequency sharp nebulous 
features in particular those varying on scales less than the 
median background filter. On this, the standard photomet¬ 
ric aperture (of 1"), and the associated background annulus 
of each source were marked. Any sources with non-uniform 
filamentary structures or nebulosity within the background 
annulus were excluded from further study. It is noted that 
previous similar work within the Lagoon Nebula has faced 
similar difficulties, including the slit spectroscopic work of 
Arias et al. (2007) where the authors were unable to subtract 
accurately the background emission in select cases, and in 
the fibre spectroscopic work of Prisinzano et al. (2007; 2012) 
where the authors were unable to measure accurately Ha 
equivalent widths of stellar sources due to difficulty in isolat¬ 
ing the background nebulosity from stellar emission. In par¬ 
ticular, the bright and complex nebulosity around the Hour¬ 


glass nebula (a 2. &' x 2.6 / region centred on RA 18 h 03 m 41 s , 
Dec —24°22 , 18") varies on scales less than a few arcsecs, and 
makes it extremely difficult to isolate point-source Ha emis¬ 
sion. Hence, all Ha sources within this region were excluded. 
Based on the single-slit observations of Arias et al. (2007), 
we estimate that up to eight known accreting PMS visible 
in the optical may have been omitted because of removing 
this region. 

In the resulting sample, 12,384 sources have riHa pho¬ 
tometry. A subset of 11,060 have gr-band photometry and 
4,864 have u-band photometry. The selected sources are 
shown in the (g — i ) versus i CMD in Fig. [2] The saturation 
limit of the VPHAS+ photometry (g ~ 14 mag) means that 
an upper main sequence containing the bright OB stars of 
the Lagoon are saturated is not visible. The lack of faint 
background stars first remarked by Prisinzano et al. (2005) 
is also seen, and is interpreted as a blue envelope due to 
the molecular cloud preventing the detection of background 
stars at visible magnitudes (see Section 1). We also note 
a sequence of reddened late-type dwarf stars in the range 
(3.2 < (g — i) < 5) and (17 < g < 22) also found in the CMD 
of Sung et al. (2000) and Prisinzano et al. (2005). 


3 IDENTIFYING CLASSICAL T TAURI 

CANDIDATES USING PHOTOMETRIC Ha 
EQUIVALENT WIDTH 


The (r-Ha) colour is a measure of Ha line strength relative 
to the r-band photospheric continuum. Main-sequence stars 
do not have Ha in emission. Modelling their (r— Ha) colour 
at each spectral type allows for a template against which 
any colour excess due to Ha emission can be measured from 
the observed (r-Ha) colour. The (r-Ha) colour excess (de¬ 
fined as (r-Ha)excess = (r-Ha)observed-(r-Ha)modei) can 
be used to compute the Ha equivalent width (EWhc*) (see 
De Marchi et al. 20101. 


We use the (r — i) colour as a proxy for the spectral 
type. This is a reasonable assumption for late spectral types 
(K-M6) which form the majority of our observed sample 
given the 1.5 mag spread in the (r — i) colour across this 
spectral range, and a mean (r — i ) colour error due to the 
combination of random photometric noise and extinction un¬ 
certainty < 0.15 mag. This results in the average error on the 
estimated spectral type caused by using the (r — i) colour as 
a proxy for spectral classification to be between 2-3 spectral 
subclasses. To calculate the main-sequence colours, |Pickles| 
(1998) spectra for 05-M6 spectral types were convolved with 
the appropriate filter band passes and CCD response. 


In Fig. [3] the (r— Ha) vs. (r — i) diagram is plotted. The 
dashed line represents the interpolated model colours red¬ 
dened by the mean reddening E(B — V) = 0.35 (see Sec. 1). 
At a given (r—i), the (r— Ha) eX cess for each star is computed. 
The EWhc is given by: 


EWho = W x [1 - l 0 °- 4x(r — Ha) “ oeM ] (1) 


as shown by De Marchi et al. (2010). W is the rectangular 
bandwidth of the Ha filter. The photometric EWho for all 
stars having r*Ha magnitudes in our sample are measured. 
We validate our measured EWhh against spectroscopic mea¬ 
surements in Section 3.1. 
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(r - z) 



Figure 3. (i—Ha) vs. (r — z) diagram. The dashed line is the in¬ 
terpolated model track reddened by Ay = 1. Grey dots are stars 
meeting our selection criterion, and blue crosses are stars meet¬ 
ing our selection criterion that are identified as candidate CTTS 
based on their EWe a . A reddening vector for Ay = 1 is shown. 


We homogeneously select CTTS from the sample of 
stars having measured EWh q using the spectral type-EWHa 
criteria of White & Basri (2003). We note that the median 
EWho errors due to photometric and reddening uncertain¬ 
ties vary between ~ 5 - 7 A depending on spectral type, so 
we empirically adjust the selection criteria of White & Basri 
(2003) by the same amount. We consider CTTS as those 
having spectral type F5-K5 and EWh q <- 10 A, K5-K7.5 
and EWhc* < —12 A and M2.5-M6 and EWh q < -25 A. 
The results of our selection procedure are shown in Fig. [4] 
The selection criteria used are designed to exclude chro- 
mospherically active late type stars, whose maximum mea¬ 
sured EWho ranges from —11A for a M2V spectral type to 


-24 A for an M6V spectral type ( |Barrado y Navascues fc 
Martin 20031. We select 235 candidate CTTS on this basis, 


which constitute the sample discussed in the rest of this pa¬ 
per. The photometric properties of the selected CTTS are 
given in Table 2. 


3.1 Comparison of photometric and spectroscopic 
EW Ha 

To test the accuracy of our photometrically determined 
EWhq, we compare with spectroscopically measured EWhq 


from Arias et al. (2007). Arias et al. (2007) obtained inter¬ 
mediate resolution spectra covering the Ha line from the 
Magellan telescope (Las Campanas Observatory, Chile) of 
46 stars in the Lagoon Nebula. Stars were selected for ob¬ 
servation on the basis of either Ha excess in RI Ha pho¬ 
tometry from Sung et al. (2000), disc excess in near-infrared 
JHKs photometry, or extended or knotty morphological ap¬ 
pearance in Hubble images. Out of the 46 stars observed by 
Arias et al. (2007), 37 were identified as Ha emission line 
stars. Twenty seven stars have VPHAS+ rzHa photometry. 
None of the ten stars classified as Ha emission-line stars that 
are not detected in VPHAS+ Ha photometry have literature 
Ha photometry in ground-based imaging from Sung et al. 
(2000). Of the ten, six were selected on the basis of their 
morphological appearance in Hubble imaging by Arias et al. 


Figure 4. The results of the spectral type-EWna selection cri¬ 
teria. Ordinate shows the (r — i) colour with the corresponding 
spectral type labelled, while the abscissa shows the EWho of se¬ 
lected CTTS. The median error corresponding to a given EWh q 
is shown at (r — z) = 0.45. The solid line is the empirical selection 
criterion of White & Basri (2003), and the red dashed line is the 
spectral type-EWHa criteria of Barrado y Navascues and Martin 
(2003) is shown for comparison. 



Figure 5. Comparison of the photometrically determined EWh 0 
(ordinate) and spectroscopically measured EWh 0 (abscissa) from 
Arias et al. (2007) for the stars selected as CTTS (crosses). The 
solid line represents the 1:1 relation, and the dotted lines represent 
the 2a error bars. Dashed lines represent EWn a = 10 A 


(2007) and cannot be resolved to point sources in VPHAS+ 
Ha imaging (ID. 3-8 inclusive in Table 2 and 3 of Arias et al. 
2007; all located in the Hourglass nebula), or ground-based 
Ha imaging in Sung et al. (2000). Two are surrounded by 
intense nebulosity and cannot be resolved as point sources. 
Two sources were selected for observations by | Arias et al.| 
120071 using near-infrared photometry and do not have re¬ 
solvable optical counterparts (ID. 1, 9, 10, 11, in Table 2 and 
3 of Arias et al. 2007; all located in the Hourglass Nebula). 
A further two stars having rzHa photometry are saturated, 
and their photometric EWho are possibly incorrect (ID. 15, 
45 in Table 2-3 of Arias et al. 2007). 


We list the spectroscopic and photometric EWho of the 
remaining stars in Table 3. We stress that this comparison 
includes all stars having rzHa photometry for completeness. 
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Table 2. VPHAS+ photometry of candidate CTTS in the Lagoon Nebula. The full table is available in the online version of this journal. 


JNAME 

RA 

(J2000) 

Dec 

(J2000) 

r 

(mag) 

(r - i) 

(r—Ha) 

(9 - r ) 

O ~ 9) 

18023206-2415320 

270.63358 

-24.25888 

14.83i0.001 

0.63i0.002 

0.45i0.002 



18023355-2414548 

270.63977 

-24.24856 

16.40i0.003 

1.OliO.004 

0.67i0.006 



18023368-2418022 

270.64032 

-24.30061 

19.96i0.049 

1.95i0.051 

1.98i0.067 



18023560-2413024 

270.64832 

-24.21733 

18.29i0.012 

1.44i0.014 

1.17i0.023 

1.89i0.029 

0.16i0.065 

18023730-2416242 

270.65543 

-24.27338 

16.51i0.003 

1.13i0.004 

0.81i0.006 



18023836-2419313 

270.65985 

-24.32536 

18.52i0.014 

1.13i0.018 

0.88i0.032 

1.03i0.019 

—0.21i0.026 

18023897-2414274 

270.66238 

-24.24094 

17.81i0.009 

0.90i0.012 

1.24i0.014 



18023972-2419310 

270.6655 

-24.32527 

18.45i0.014 

1.31i0.016 

0.92i0.027 

1.76i0.036 

0.40i0.104 

18024000-2419346 

270.66666 

-24.32627 

17.67i0.007 

1.16i0.009 

0.85i0.014 

1.63i0.015 

0.36i0.039 

18024097-2412163 

270.67072 

-24.20452 

16.59i0.004 

0.96i0.005 

0.82i0.009 

1.43i0.006 

0.25i0.013 


Table 3. Comparison of detected Ha emission line stars with spectroscopic EWh q from Arias et al. (2007; A07), and archive data from 
Sung et al. (2000; S00). 


ID 

r 

r —Ha 

V 

R—Ha 

Sp. EW Ha 

EW Ha 

CTTS 

Comment 

A07/S00 

(mag) 


(mag) 


(A) 

(A) 



12/174 

16.27 

1.05 

16.90 


— 15.1 i 0.1 

-73.5 i 12 

N 

No lit. Ha, near HG 

Does not meet photometric criteria 

13/184 1 

15.66 

1.44 

16.53 

-2.88 

— 126 i 9 

-138.3 i 30 

Y 


15/240 1 

10.63 

0.54 

11.72 

-4.23 

-24.1 i 0.4 

-46.8 i 2 

N 

Saturated, HAeBe star 

17/292 1 

16.25 

0.95 

16.84 

-4.06 

-52 i 3 

-58.3 i 11 

Y 


18/388 

15.07 

0.57 

16.70 


-3.8 i 0.2 

— 12.8 i 2 

Y 

No lit. Ha 

19/390 1 

15.22 

0.54 

15.58 

-3.75 

—16 i 0.6 

-33.5 i 5 

Y 


21/418 1 

14.81 

0.6 

16.5 

-3.98 

-5.6 i 0.03 

— 13.4 i 1 

Y 


22/422 

14.78 

0.82 

15.53 


-27.8 i 0.5 

-35.8 i 3 

Y 

No lit. Ha 

24/440 1 

14.1 

0.44 

15.99 

-4.14 

-0.65 i 0.03 

-0.75 i 1 

N 

Does not meet EWn Q -Sp. type criteria 

25 2 

16.47 

0.46 



-8.1 i 0.3 


N 

No lit. Ha, no Ha emission detected 

26/482 

13.95 

0.50 

15.45 

-4.22 

-9.3 i 0.2 

-8.6 i 1 

N 

Does not meet EWn Q -Sp. type criteria 

27/486 1 

16.79 

1.14 

16.9 

-3.26 

-75 i 4 

—82 i 17 

Y 


29/495 1 

15.89 

0.72 

16.73 

-3.80 

-29 i 1 

-22.5 i 4 

Y 


30/508 1 

16.44 

1.54 

16.73 

-3.73 

—30.3 i 0.7 

— 155.7 i 33 

Y 

Marked as eruptive variable star 
in Kukarkin & Kholopov (1982). 

31/531 

16.29 

0.59 

16.03 

-4.18 

-14.96 i 0.04 

— 13.10 i 3 

Y 


32/540 1 

15.56 

0.61 

16.17 

-4.16 

—6.7 i 0.4 

-10.1 i 2.5 

Y 


33/547 1 

14.33 

0.6 

14.93 

-4.16 

-2.26 i 0.07 

-17.4 i 2 

Y 


34/552 1 

14.55 

0.71 

15.39 

-3.80 

— 12.8 i 0.3 

—49.7 i 7 

Y 


35/666 1 

15.68 

0.68 

16.53 

-3.90 

—35.4i 0.3 

-20.5 i 2 

Y 


36/681 

13.37 

0.39 

15.0 

-4.59 

— 0.7 i 0.07 


N 

No Ha emission detected 

40/726 1 

16.07 

0.8 

16.88 

-3.91 

-48 i 2 

-45.9 i 3 

Y 


41 2 

13.20 

0.35 



-3.18 i 0.06 

-2.09 i 1 

N 

No lit. Ha 

Does not meet EWnct-Sp. type criteria 

42/782 1 

15.47 

0.87 

16.25 

-3.55 

—54.0 i 2 

-50 i 4 

Y 


43/810 1 

16.29 

2 

16.25 

-3.96 

—20.8i 0.3 

-22.3 i 2 

Y 


44/862 1 

15.34 

0.96 

16.90 

-3.26 

— 22.3 i 0.4 

—41.8 i 4 

N 

Does not meet photometric criteria 
Located in M8 E rim 

Outlier in EWjja v s. (R— Ha) 
diagram in Arias et al. (2007; Fig. 5) 

45/879 1 

10.74 

0.93 

11.94 

-3.15 

-80 i 1 

— 101.8 i 10 

N 

Saturated, HAeBe star 

46 2 

15.07 

0.54 



-2.85 i 0.2 

— 19.1 i 3 

N 

No lit. Ha 

Does not meet photometric criteria 


1 Object selected as PMS stars with strong Ha emission in Sung et al. (2000); 2 Object not detected in Sung et al. (2000) 


Sources that are not selected as CTTS either because they 
do not meet the EWhq -spectral type selection criteria, or 
photometric quality checks are marked. The sample of se¬ 
lected CTTS does not include stars that have not met these 
criteria. We plot the spectroscopic and photometric EWhq 
of the stars meeting the quality criteria in Fig. [5] We End a 


good correlation between the spectroscopic and photomet¬ 
ric EWho for all objects with a median scatter ~8A, ap¬ 
proximately the error on the photometric EWho due to the 
assumption of uniform extinction and random photometric 
noise. The errors dominate the scatter for objects meeting 
the quality criteria and with EWho > —10 A, which is the 








Ha and u-band accretion rates in the Lagoon Nebula 7 



Figure 6. Unsharp mask Ha thumbnails (3" X 3") of the three 
CTTS (bottom panel, upper right panel; IDs. A0719,30,34 in 
Table 3), and a Ha emission line source that is affected by variable 
nebulosity (upper left panel; IDs A07 44 in Table 3) all which 
have photometric EWh q 2 a greater than spectroscopic EWh q . 
The solid line marks the photometric aperture, and the dashed 
line the background correction aperture. 



r (mag) 


Figure 7. r-band histogram of CTTS candidates identified in 
this work (blue bars) overlaid by Sung et al. (2000) CTTS can¬ 
didates (grey bars). The majority of new CTTS are identified at 
r < 17 mag, beyond the limiting magnitude of Sung et al. (2000), 
while a significant fraction at magnitudes smaller than this. This 
is because our EWh 0 identification method can separate CTTS 
lying near the main-sequence, which do not have large Ha ex¬ 
cesses. 


minimum threshold used in this study, because small pho¬ 
tometric or reddening uncertainties will produce relatively 
large errors in EWho for stars with weak-Ha emission, as 
they lie close to the main sequence. Three CTTS show dis¬ 
crepant EWho (> 2<t). It has been observed that the EWh q 
of a small fraction of accretors in a given sample can vary 
between —7 to —100 A on time-scales greater than a year 
(Costigan et al. 2012) due to natural variations in the accre¬ 
tion process. This may be apt in the case of A07 30/ S00 508 
(see Table 2) as it has been marked as an eruptive vari¬ 
able (GCVS V1780 Sgr) in |Kukarfri n & K holopov] (1 982| . 
Therefore, we suggest that the difference between the spec¬ 
troscopic and photometric EWhq is most likely due to in¬ 
trinsic accretion variability. To demonstrate that the pho¬ 
tometric quality checks ensure that nebulous emission does 
not affect the measured EWhq of the selected CTTS sig¬ 
nificantly, we checked the Ha images of stars whose spec¬ 
troscopic and photometric EWho differ ~ 2cr, showing their 
object and background annuli overlaid (see Section 2.1) to 
ascertain the absence of rapidly spatially varying nebulosity 
in the background annulus (Fig. [d| on scales smaller than 
can be corrected for. For comparison a star which is affected 
by such nebulosity is also displayed. Finally, although steps 
have been taken to ensure that the results are not affected 
by nebulosity, we accept that due to the difficulty of separat¬ 
ing the strong sky nebulosity from stellar emission especially 
in Ha photometry, some residual may be present as is also 
the case in spectroscopic work in the region (Arias et al. 
2007; Prisinzano et al. 2007; 2012). Due to the small num¬ 
ber statistics of comparable spectroscopic EWho, it is not 
possible to make a statistical statement. We therefore wel¬ 
come further comparisons of the photometric EWho with 
spectroscopic measurements, particularly from future large 
surveys such as the GAIA-ESO Survey (GES; Prisinzano et 
al. in prep.) 


3.2 Comparison with literature Ha emission line 
star lists 


We estimate the recovery rate of CTTS candidates identi¬ 
fied using our method, and thereby the validity of our sam¬ 
ple for a statistical analysis by comparing with all available 
literature optical studies of Ha emission line stars in the re¬ 


gion (Sung et al. 2000 Prisinzano et al. 2007 Arias et al. 


|2007[ Prisinzano et al 7| |2 012 [ ). |Sung et a.l.| ( |2000[ ) observed 

the central NGC 6530 region (a ~20' square centred on RA 
18'*04 m 09 3 , Dec -24°21'03") in UBVRMa in a study lim¬ 
ited to V < 17 mag. They identified 37 stars from the (R— Ha 
)vs. (V — I) diagram as stars having (R— Ha) excess > 0.36, 
which they classified as accreting PMS stars (CTTS). 

All 37 stars identified as CTTS candidates by Sung et al. 
(2000) have VPHAS+ photometry, although four sources are 
saturated and their photometric EWho may be inaccurate. 
The remaining stars have photometric EWh« <— 8 A indi¬ 
cating that they are reliably identified as emission-line stars 
using VPHAS+ photometry. Additionally, in the region of 


the sky observed by Sung et al. (20001 we identify 125 new 


CTTS candidates. This is roughly three times as many as 
identified by Sung et al. (2000). A comparison of the magni¬ 
tude distribution of the two samples is plotted in Fig. [T] We 
converted Sung et al. (2000) RI magnitudes to VPHAS+ 
r-band using the colour equations of jJester et a l.] ( |2005| ). A 
significant fraction of VPHAS+ identified CTTS candidates 
are located at r > 17 mag, beyond the limiting magnitude of 


Sung et al. (2000). We roughly double the number of Sung 


et al. (2000) CTTS candidates at r < 16 mag, as we can sep¬ 
arate them from main sequence stars on the basis of their 
EWho, without needing to rely on large (R— Ha) excesses. 

Prisinzano et al. ( |2007[ ) presented Ha fibre spectra of 
95 stars in the Lagoon Nebula, but experienced difficulty in 
separating the stellar Ha component from any nebular con¬ 
tamination. They do not measure stellar EWho to classify 
accreting PMS stars, instead employing the Full Width Zero 
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Ha FWZI Prisinzano et al. 2007; 2012 (A) 


Figure 8. FWZI of Ha line profiles measured by Prisinzano et 
al. (2007; 2012) plotted against EWn n measured using the pho¬ 
tometric method. All of the stars plotted were classified as CTTS 
by Prisinzano et al. (2007;2012) based on the measured Ha FWZI 
and line profile. We find that almost all these stars are in emission 
in the VPHAS+ sample, but a significant majority of them are 
below the minimum EWh q threshold for CTTS (dashed line). 
We find no relation between Ha FWZI and EWh q - 


Intensity (FWZI) of the Ha line to do so. They argue that 
broad emission due to accretion can be separated from any 
narrow nebular contamination using the FWZI parameter. 


Prisinzano et al. (20071 classify 31 of the 95 stars for which 
they present Ha spectra as CTTS (Fig. [8} purely based on 
the FWZI value if the Ha line is in emission. 30 are identified 
as having Ha in emission in our sample, although slightly 
less than half are below the spectral type-EWHa selection 
criteria we apply to identify CTTS. We find no correlation 
between the Ha FWZI and EWhc«. Prisinzano et al. (20121 


used new Ha multi-slit observations to add a further 16 
accreting stars to the sample of |Prisinzano et al. (20071, 
which are also shown in Fig. [8] All of these stars have Ha 
in emission, but around a third of the sample fall below the 
spectral type-EWim selection criteria. In short, our method 
can reliably identify Ha emission in all stars that have been 
previously classified as Ha emission-line stars using spectra. 
We stress that although the FWZI of the Ha line is a good 
identifier of accretion, but may not be an accurate measure 
of the accretion rate in CTTS ( White fc Basri||2003 Cost! 
gan et al.|2014 l. 


The comparison of our sample with literature Ha emis¬ 
sion line star lists suggest that all known CTTS sources 
within the VPHAS+ dynamic range can be reliably re¬ 
covered. We identify a total of 235 CTTS candidates over 
our area of study, which contains at least 170 new accret¬ 
ing sources that represent ~ 200 per cent increase in known 
CTTS stars within the region. The area covered by this 
study is significantly larger (around 3.5 times) than previous 
photometric searches for PMS stars in the Lagoon Nebula 
(see Tothill et al. 2008 for a summary). We conclude that the 
uniform expansive sky coverage, and photometric depth af¬ 
forded by the VPHAS+ survey can homogeneously identify 
a substantial number of CTTS candidates; and the photo¬ 
metric method estimates the EWho accurately in compari¬ 
son with intermediate-resolution spectroscopy. 


4 RESULTS 


In the current accretion paradigm, it is assumed that the 
stellar magnetosphere truncates the circumstellar disc at an 


inner radius (Hartmann 1994 Muzerolle, Calvet & Hart- 
|mann 1998}). Gas flows along the magnetic field lines at this 
truncation radius, releasing energy when impacting on the 
star. It is generally accepted that the accretion energy goes 
towards heating and ionising the gas, and that the accretion 
luminosity (L acc ) can be measured from the reradiated emis¬ 


sion line luminosity (Calvet & Gullbring 1998 Hartmann 
2008). The M acc can then be estimated from the free-fall 


equation: 


M acc = 


LaccR* 

CM, 


Rii 


Rin — R* 


( 2 ) 


M* and R * are the stellar mass and radius respectively. Ri n 
is the truncation radius. 

In this Section, we estimate the stellar and accretion 
properties of our candidate CTTS based on the magneto¬ 
sphere accretion paradigm. In Section 4.1 we estimate their 
stellar masses and ages. We estimate their M acc based on 
their EWho in Section 4.2 and it-band excess in Section 4.3. 
We discuss their near-infrared properties in Section 4.4. The 
caveats, completeness and contamination of our sample are 
discussed in Sections 4.5-4.8. 


4.1 Stellar properties 


4-1.1 Mass and age 

The mass (A/,) and age (i„) of each star is estimated by 
interpolating its position in the observed r versus (r — i) 
CMD with respect to PMS tracks and isochrones (Fig. [9}. 
Siess, Dufour & Forestini (2000) solar metallicity single star 


isochrones and tracks were employed as they cover the mass 
and age range of our sample. The isochrones were reddened 
assuming the mean reddening E(B — V ) = 0.35 following a 
reddening law Rv =3.1 (see Section 1). The systematic dif¬ 
ferences that would arise on adopting different stellar mod¬ 
els are discussed in Section 4.6.1. The tracks and isochrones 
were converted from the luminosity- temperature plane to 
the colour-magnitude plane using the colour equations of 


Jester et al. (20051. 


Fig. 10 a) shows the luminosity function of our sample. 


The interpolated mass distribution of our candidate CTTS 
is plotted in Fig.|10[b). All our candidate CTTS have masses 
between 0.2-2.2 Mq. The mass distribution shows the ex¬ 
pected increase towards lower masses, peaking at 0.35 M. 
The mass cut-off and steady decline in sub-solar mass below 
0.3 Mq is expected because of the limiting magnitude. 

The ages of stars in our sample range between 0-9 Myr, 
with 95 per cent of stars having an age less than 2 Myr 
(Fig. [To] b). Few stars fall beyond the 30000 yr isochrone. 
For these stars, age estimates represent an upper limit. The 
median age of our sample is 0.85 ± 0.6 Myr. This agrees with 
the main-sequence lifetime of the early-type stars found in 
the region (Tothill et al. 2008}. In particular, 9Sgr (04V; 
Walborn 19731 and H36 (07V; van Altena & Jones 1972) 
are two prominent O-type members still lying on the zero- 
age main sequence (ZAMS), and are expected to evolve off 
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the ZAMS after 1-3 Myr. Chen et al. (2007) measured an 
intrinsic velocity dispersion of 8 kms“ . Assuming this as 
an expansion velocity, Tothill (2008) calculated a dynamical 
age of 1 Myr. Furthermore, the age of the PMS stars that 
have X-ray emission (Damiani et al. 2004) or Ha emission 
(Sung et al. 2000; Arias et al. 2007) are estimated to be 
around 1 Myr. 

We caution the reader that higher mean age estimates 
(> 2 Myr) for the general population of stars in the Lagoon 
Nebula have been measured in the literature. In particular, 
Walker (1957) and van Alten & Jones (1972) suggest an age 
estimate ~ 3 Myr, although both of these studies only fit 
a ZAMS and use a distance modulus of around 11.5 mag. 
Recent works by Mayne et al. (2007) and Mayne & Nay¬ 
lor (2008) use the CMD diagram to build a relative age 
ladder, and obtain an age of 2 Myr. Also, van den Ancker 
et al. (1997) find a highly probable member with an age of 
15 Myr and suggest a prolonged period of star formation tak¬ 
ing place within the Lagoon Nebula. The evidence for this 
remains unclear. The study by Prisinzano et al. (2005) dis¬ 
cussed the possibility of an age spread, by comparing their 
optical sample with the X-ray identified PMS stars of Dami¬ 
ani et al. (2004). They find a median age of 2 Myr using their 
optical sample, with the older stars preferentially located to 
the northeast, and younger stars located to the south and 
west. In our sample, higher or lower ages may be estimated 
for a few stars because the assumption of uniform distance or 
extinction may be inappropriate, or because they are non¬ 
members and not representative of the age of the sample. 
Since our study is limited to an assumption of uniform ex¬ 
tinction and no strict membership tests barring the EWhq 
criteria, we cannot differentiate between ages that are in¬ 
correct due to these effects or genuinely older. Therefore, 
we cannot confirm or rule out the presence of a true age 
spread, and further data on cluster membership is required 
to explore this possibility. 

Finally, we stress that the precise ages of individual 
stars derived using isochrones may be uncertain, as most of 
our stars are around one million years old. At this point in a 
star’s evolution, age determinations using stellar isochrones 
are imprecise because the age is thought to be still depen¬ 
dant on the accretion history and stellar birthline correc¬ 
tions. However, the median age of the CTTS population 
(~ 1 Myr) is relatively accurate given large number statis¬ 
tics, and the fact our result agrees well with age estima¬ 
tions based on the expansion velocity (Chen et al. 2007) 
main-sequence lifetime of OB stars, and literature age es¬ 


timations of the pre-main sequence population (Sung et al. 
2000 Damiani et al.|[2004 Arias et al.|2007 l. 


4.2 Accretion properties 

4-2.1 Ha line luminosity 

The Ha line flux (Fu a ) is calculated by subtracting the stel¬ 
lar continuum flux (-F corl tiimum) from the total flux in the 
Ha band (Ftotai). The total Ha flux (F t ota i) is given by the 
unreddened Ha magnitude (rriHa) according to the relation, 


E tota l = F 0 X [ 1Q —°- 4(mH “ + 0.03)]. 


(3) 


m-Hc, is the unreddened Ha magnitude. Fq is the band- 
integrated reference flux. For the VPHAS+ Ha filter, 


Fo=l.84x10 7 ergscm 2 s 1 . Ftotai is a sum of the contin¬ 
uum (^continuum) and line flux (Fka)- 

File — F total Fcontinuum* (4) 

Following from equation 1, equation 4 can be written as 

-EWho/W 


Fua — Ttotal X n ,,, /xxt ' 

1 - EW Ha /W 

The line luminosity is given by Lhu'- 

Ltso, = Fhc x 4-7T d 2 . 


(5) 


( 6 ) 


Here the distance d = 1250 ± 50 pc. Contamination from 
N[n] AA 6548, 6584 is excluded by assuming it causes 2.4 
percent of the Ha intensity (De Marchi et al. 20101. 


4-2.2 Accretion luminosity (T/ ac< J 

Measurements of L acc from UV continuum spectral mod¬ 


elling in the literature (Sicilia-Aguilar, Henning & Hartmann 


2010| |Herczeg fe Hillenbrand|2008| |Dahm|2008 Hartigan & 


Kenyon|2003| ) were plotted against Ln a by Barentsen et al. 

(2011 1 . A power law relationship, L acc ex LuH was found, 


having an index 7 = 1.13, consistent with the accretion mod¬ 
els of |Muzero llc et al. (i l998| . Based on results of Barentsen 
et al. ( 2011 ), we assume: 

log L acc = (1.13 ± 0.07)log L Hq + (1.93 ± 0.23) (7) 

log L acc and log Lu a are measured in units of solar lumi¬ 
nosity, Lq. The root mean scatter = 0.54. The scatter in 
the observed relationship is likely caused by a combination 
of circumstellar absorption, uncertain extinction and line 
emission caused by processes other than accretion. 


4-2.3 Mass accretion rate A7 ai 


Ms.cc are estimated from Equation 2. We adopt 
Rixe = 5 ± 2 from ( |Gullbring et ah] |1998| |Vink et al.| 


2005). The res ultan t distribution of mass accretion rates are 
plotted in Fig. 10 d). The median M acc is 10 -8 ' 1 M^yM 1 . 


Our results are given in Table 4. 


4.3 Comparison between u-band and Ha mass 
accretion rates 

91 candidate CTTS have ugr photometry from the blue data 
set meeting our selection criteria. Four stars do not show any 
measurable UV-ex cess in the (u — g) vs. (g — r) diagram, 
and have EWm* > —20 A. We discuss them in more detail in 
Section 4.8. For the remaining stars, we measure M acc from 
the excess ultraviolet emission (Gullbring et al. 19981. 

Similar to the method described in Section 3, the (g — 
r) colour is used to estimate spectral type, while the (u — 
g) colour is the measure of UV-excess due to the accretion 


process. The UV-excess stars are shown in Fig. 11 a). The 
excess flux in the u-band (F UtexC e S s) is: 


- u, excess 


= F 0 , v X [10~ U ° /2 ' 5 - 1 Q-(—9)»od.l+90)/2. 5 ]^ (g) 


Fo, u is the u-band integrated reference flux, u 0 and go are 
the observed dereddened magnitudes, and (u — g) m odei is the 
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Figure 9. r versus ( r-i ) CMD of CTTS candidates. Labelled isochrones (solid lines) and tracks (dashed lines) are from Siess et al. (2000). 
Model isochrones and tracks are reddened by E(B — V) = 0.35 following a reddening law Ry = 3.1. Circles indicate the position of each 
candidate CTTS, and symbol size is representative of logM acc as shown by the legend in the top-right hand corner. The reddening vector 
for Ay = 1 is shown. The median magnitude errors in a given 2 mag bin, and the associated colour errors are shown at the corresponding 
magnitude and (r — i) =0.02. 


corresponding model colour from Castelli & Kurucz (20041. 
The excess flux is used to estimate the L acc using the em¬ 
pirical relation (Gullbring et al.||1998l; 


log L acc — log Lu.excess T 0.98. (9) 

Lu,excess is the u-band excess luminosity. L acc and L„,excess 
are in units of solar luminosity L@. Mass accretion rates 
are calculated using Eq. 2, with results given alongside Ha 
derived M acc in Table 4 . T he (u — g) e xcess are compared 
with the Haexcess in Fig. |ll[ (b), and M acc thus determined 
are compared with the ones determined from the Ha-line 
luminosity in Fig. |12| 

The (u — g)excess and Ha e xcess, two independent (ex¬ 
cept for the assumed extinction value) measures of accre¬ 
tion correlate well. We find the best fit relation to follow 

Ha 

excess — 0.3 (u — g) excess 0.14 (Fig. |ll| j). There are only 
six low accretion rate objects (logM acc < —8.5) detected in 
the u-band, that are detected in Ha (Fig. |10| i). This is due 
to the faint magnitude limit in u ~ 19 mag, thereby exclud¬ 
ing the low-mass stars detected in Ha. The logarithm of 
the mass accretion rates determined from the Ha-line lu¬ 
minosity and u-band excess luminosity have a mean varia¬ 
tion (R.M.S. scatter) of 0.17dex. There exist a minority of 


CTTS candidates (four) which have differences in logM acc 
~0.8dex. Since the observations are non-contemporaneous, 
the M acc should display intrinsic variabilities over long time- 
scales with maximum variations in log M acc expected to be 
<0.52dex (Cost igan et al.|2012 |. 


This result indicates no significant statistical effect of 
accretion variability on the sample as a whole. However, 
accurate u band photometry of low-mass stars is compara¬ 
tively more challenging to obtain. Therefore Ha photometry 
is better suited to provide statistically significant samples of 


4.4 Infrared properties 


A simple sanity check on the validity of the identified CTTS 
can be made by studying their near-infrared properties. This 
is because CTTS display near-infrared excesses when com¬ 
pared to a purely stellar template due to the presence of dust 
in their inner circumstellar discs (Cohen & Kuhi]l979). 

Near-infrared JHK photometry in the Lagoon Neb¬ 
ula from the UKIRT Infrared Deep Sky Survey (UKIDSS) 
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Figure 10. Distribution of magnitude M r (a), Mass (b), Age (c) and M acc (d) of our sample. 




Figure 11. (a) ( u — g) versus ( g — r ) diagram. The dashed line is the reddened model colour track from Castelli &; Kurucz (2004). Stars 
meeting our selection criteria (grey dots) and CTTS identified on the basis of their EWh 0 (blue crosses) are shown. A reddening vector 
for Ay = 1 is shown, (b) Comparison of ( u — g) e x cess and Ha excess of stars in (a). The solid line represent the best-fit relation having 
HcKexcess — 0.3 ('ll . 7 )excess 0.14. 


Galactic Plane survey (Lucas et al. 20081 was cross-matched 
with the positions of the CTTS in our sample. Only cross¬ 
matches having random error <0.1 mag in JHK, and clas¬ 
sified as having a stellar profile in the UKIDSS catalogue 
were selected. UKIDSS reaches a depth of 5cr at K ~ 19 mag, 
which translates to roughly a star of spectral type M6 in 
the Lagoon Nebula, suggesting that most of the stars in our 
sample should be identifiable in the UKIDSS survey. 224 
(95 per cent) stars in our sample were found to have UKIDSS 
cross-matches. 

In Fig. |13[ a), the (J — H) vs. (H — K) colour-colour 
diagram of stars located in the Lagoon nebula is plotted. 
Also shown are near-infrared counterparts of cross-matched 
candidate CTTS. The dashed line is the CTTS locus of 


cess emission using disc accretion models having log M acc 
between —6 to —8 MQyr~ . The solid line is the main se¬ 


Meyer, Calvet & Hillenbrand (1997) which predicts the ex¬ 


quence locus from Bessell, Castelli & Plez (1998). Most stars 
(~ 85 per cent) having near-infrared colours in our sample 
lie on the CTTS locus, suggesting on the basis of the near- 
infrared diagram alone that they are CTTS. This confirms 
our sample as consisting primarily of CTTS undergoing ac¬ 
cretion. Around 35 stars lie near the main sequence locus, 
indicating that they are either weak-line T Tauri stars, or 
have no excess infrared emission suggestive of a circumstel- 
lar disc. Three of the stars lying below the main sequence 
spur are the stars that showed no UV-ex cess. We discuss 
them further in Section 4.8. 

We also cross-matched our sample with mid-infrared 
Spitzer data (Kumar & Anandarao 2010). 32 cross- 
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Table 4. Derived properties of candidate CTTS in the Lagoon Nebula. The full table is available in the online version of this journal. 


JNAME 

EWhc, 

(A) 

Mass 

(Mg) 

Log Age 
yr 

^acc,He* 
(Mgyr -1 ) 

-M a cc, U 
(Mg yr- 1 ) 

18023206-2415320 

-14.34T2.66 

l.OiO.OO 

5.6i0.08 

-7.7T0.6 


18023355-2414548 

-18.90i3.27 

0.4i0.01 

4.7i0.19 

— 7.9i0.6 


18023368-2418022 

—213.96i36.9 

0.2i0.20 

5.5i2.54 

-8.0T0.6 


18023560-2413024 

—63.16il0.8 

0.3i0.05 

5.4i0.69 

—8.1i0.6 

—7.8T0.7 

18023730-2416242 

— 28.39i5.15 

0.4i0.01 

4.6i0.20 

— 7.6i0.6 


18023836-2419313 

—36.71i4.95 

0.4i0.06 

6.2i0.85 

—8.8i0.6 

-7.7T0.7 

18023897-2414274 

— 110.99i21.0 

0.5i0.04 

6.3i0.54 

—8.1i0.6 


18023972-2419310 

—33.25i4.58 

0.3i0.06 

5.6i0.77 

—8.6i0.6 

—8.2T0.7 

18024000-2419346 

—31.36i5.14 

0.4i0.03 

4.7i0.43 

—8.3i0.6 

-7.9T0.7 

18024097-2412163 

—39.45i7.23 

0.5i0.01 

5.7i0.23 

— 7.7T0.6 

— 7.2T0.7 



Figure 12. Comparison of M ac c derived from Ha line luminosity 
and u-band excess luminosity. Mean error bars are shown in the 
upper left hand corner. 


matches were found. In Fig. [ML)’ the 3.6 fi.m—4.5 pm vs. 
5.8 |im—8 pm colour-colour diagram is plotted. The solid 
box represents the typical colours for PMS sta rs having 
mean M acc = 10 -8 Mgyr -1 ( D’Alessio et al.||200l K Most of 
the cross-matched CTTS candidates in our sample lie in 
this box, providing an additional sanity check on the re¬ 
sults. Six stars lie at a 5.8 pm-8 pm colour much brighter 
than maximum expected colour. Their accretion rates are 
approximately a magnitude higher than the model at 
~ 10“ 7 M 0 t/r _1 . 


Table 5. Errors in the derived parameters. 


Parameter 

Median 

Range 

<j 

ew Ho (A) 

-47.0 

-11 to -234 

7.0 

Mass (Mq) 

0.4 

0.2 to 2.2 

0.07 

Age (Myr) 

0.85 

<0.1 to 9 

0.6 

logM acc (Mgj/r^ 1 ) 

-8.1 

—6.3 to —9.4 

0.6 


The true distance for some stars may vary from the assumed 
distance. To account for the error due to the uncertainty in 
the distance, we assume a ad = 50 pc. 50 pc is the width of 
the region at the assumed distance. The RMS scatter of 0.54 
in the I/HaMacc relation was used to calculate errors aris¬ 
ing due to the deviation from this relationship for individual 
objects. 

The mean errors (a) determined are listed in Table 5. 
We find that our EWh q is accurate within 7 A while taking 
into account random photometric and extinction uncertain¬ 
ties. The stellar mass and age are constrained to 0.08 Mg 
and 0.6 Myr respectively accounting for random photomet¬ 
ric, extinction and distance uncertainties. The error on M acc 
is dominated by the scatter in the Ln a -L^ cc relation. 


4.6 Caveats 

We discuss the differences in our results arising from adopt¬ 
ing different stellar evolutionary models, and the effect of 
setting the binary fraction to zero. 


4.5 Random errors 

The random errors affecting the determined parameters are 
(i) random photometric noise (ii) scatter in the extinctions 
of individual stars compared to the adopted uniform value 
(iii) scatter in the distances of individual stars compared 
to the adopted uniform value (iv) scatter in the LHa-T acc 
relation (Eq. 7). 

The errors due to random photometric noise are prop¬ 
agated through the results. The true extinction for some 
stars may vary from the assumed uniform extinction. To 
calculate the errors arising due to the assumption of uni¬ 
form extinction, we assume a E ^ B _ V ' ) =0.2 from Section 1. 


4-6.1 Stellar evolutionary models 

The mass and radius used to calculate the mass accretion 
rates, and the age of the CTTS are determined from inter¬ 
polating their position in the CMD relative to stellar model 
tracks and isochrones. Apart from differences due to indi¬ 
vidual uncertainties in a star’s photometry or extinction, 
the model tracks and isochrones also differ between various 
authors. This is because of differences in the input physics, 
in particular the treatment of the birthline, accretion his¬ 
tory and the deuterium burning phase. Moreover, the con¬ 
version of luminosity to magnitude based on bolometric cor¬ 
rections, and the colour-temperature conversions are highly 













Ha and u-band accretion rates in the Lagoon Nebula 13 




Figure 13. (a) (J — H) vs. (H — K) colour-colour diagram. Stars located in the line of sight of the Lagoon Nebula having UKIDSS 
photometry (grey dots), and cross-matched CTTS candidates (blue crosses) are shown. The red dashed line represents the CTTS locus 
of Meyer et al. (1997). The solid line is the main sequence locus from Bessell et al. (1998). The reddening vector for Ay = 1 is shown, 
(b) The CTTS candidates having cross matches in Spitzer photometry are shown as blue crosses. The box outlines the predicted colours 
for PMS stars having M acc ~ 10 -8 Mgi/r -1 (D’Alessio et al. 2005). 


non-linear, especially for late-type stars, leading to slight 
differences. 

To estimate the uncertainties in the measured proper¬ 
ties due to differences between model isochrones and tracks 
we compare the stellar masses, and ages derived from the 
Siess et al. (2000) stellar models with the PMS models of 


Tognelli, Prada Moroni & Degl’Innocenti (2011) and Bres- 


san et al. (20121. In Fig. 14 we plot the r vs. (r — i) CMD 


of CTTS overlaid with each of the stellar models. We also 


plot the histograms (Fig. 15 I showing the estimated distri¬ 
butions of stellar masses and ages interpolated from each 
of the models compared to the distributions estimated from 
the Siess et al. (2000) models. 

We find that the median age of our sample determined 
using Tognelli et al. (2011) isochrones is 1 Myr and Bressan 
et al. (2012) isochrones 0.9 Myr. The distribution of stellar 
ages derived using the Siess et al. (2000) and Tognelli et al. 
(2011) isochrones are similar, but there is a larger number 
of CTTS having ages > 2 Myr according to the Bressan et 
al. (2012) isochrones, i.e. evidence for an age spread. This 
result is similar to that of Sung et al. (2000), where the 
authors found age spreads ranging between 2 to 5 Myr us¬ 
ing different stellar isochrones. We also point out that the 
spacing between isochrones decreases from 0.5 mag in colour 
for ages 0.1-2 Myr to less than 0.5 mag between 2-10 Myr. 
Small differences in the conversion from the theoretical to 
the observational plane of the Hertzsprung-Russell diagram 
may lead to differences in the precision of the interpolated 
ages. The turnover of the masses determined using the Siess 
et al. (2000) tracks (M*=0.35Mg) is slightly lower than 
those of Tognelli et al. (2011) and Bressan et al. (2012) at 
M*«0.4Mg. In the distribution of Bressan et al. (2012) 
masses, there is a long tail towards the higher masses. Sung 
et al. (2000) also report differences in the estimated masses 
of Lagoon Nebula PMS stars when using different stellar 
tracks. Finally, we plot the distribution of M acc calculated 
using the stellar masses and radii estimated from the dif¬ 


logMacc determined using Bressan et al. (2012) and Tognelli 
et al. (2011) are smaller by —0.05dex compared to the Siess 
et al. (2000) median of logM acc = — 8.1 Mg yr~ x . There is 
a tail of lower M acc sources found using the Tognelli et al. 
(2011) and Bressan et al. (2012) models, absent in the dis¬ 
tribution of Siess et al. (2000) M acc . The M acc estimates are 
less sensitive to the stellar models. 

Finally, we note that a significant fraction of stars with 
stellar masses M, >0.6 Mq (or r > 16.8) in Fig. [9] have com¬ 
paratively younger ages (i*<0.8Myr) than stars having 
0.2 Mg < M* < 0.6 Mg, the bulk of whose ages lie between 
0.8-2 Myr according to the Siess et al. (2000) isochrones. 
The location of these stars is similar when compared to the 


Tognelli et al. (20111 and Bressan et al. (20121 PMS mod 


ferent stellar models (Fig. 15;). We find that our median 


els (Fig. 151. It is possible that a bulk of these stars could 
be foreground reddened late-type stars. However, all stars 
having u-band photometry (53 per cent of the total sample) 
have demonstrable u-band excesses, and 86 per cent lie near 
the CTTS locus in the near-infrared colour-colour diagram. 
Evidence of multiple accretion signatures suggests that the 
number of contaminants in this sample is low. Another 
observational effect is that in the adopted EWHa-spectral 
type selection criterion, stars spanning F5-M2 spectral types 
are selected as CTTS if they have EWhc* between 10-12 
A. More massive stars have greater luminosities at similar 
EWho leading to systematically higher Ha luminosities and 
thereby accretion rates at higher masses. But, in Section 5 
our detection limits calculated based on the EWHa-spectral 
type selection criterion and stellar parameters from a 1 Myr 
isochrone indicate that the lack of low accretion rates at 
M* > 0.6 Mq is not likely to be due to the EWHa-spectral 
type criterion. When compared to other CMD of the Lagoon 
Nebula across a similar mass range, we note a very similar 
split sequence in age at 0.6 Mq in the theoretical colour- 
magnitude diagram of Sung et al. (2000; Fig. 9a in that Pa¬ 
per), and a not as prominent effect at ~ 0.8 Mq in the CMD 
of Prisinzano et al. (2005; Fig. 7 in that paper). Isochrones 
fitted to open cluster may yield good fits to one section of 
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(r-i) 


Figure 14. The r versus (r — i ) colour-magnitude of CTTS candi¬ 
dates star overlaid with the Bressan et al. 2012 (top) and Tognelli 
et al. 2011 (bottom) stellar tracks and isochrones. For both pan¬ 
els, dashed lines from bottom up represent stellar tracks of 0.1, 
0.3, 0.6, and 0.8Mq, while solid lines from right to left represent 
PMS isochrones having age 0.01, 0.3, 1, 3, and lOMyr respec¬ 
tively. 


the CMD, but may systematically deviate from the observed 


sequence in another section in some cases (Bonatto, Bica & 
Girardi||2004 Mayne et al.||2007 ). Lastly, it is possible that 


the comparatively more massive stars may have been formed 
later than the less massive stars, and the observed differ¬ 
ence in ages is real. Further information on membership is 
required to differentiate the two effects. 


4-6.2 Bias due to unresolved binaries 

Due to the unknown dependence of accretion properties on 
binarity, it is not possible to design a correction for the 
brightening effect of unresolved/unidentified binaries. We do 
not apply any corrections due to this effect. 

The angular resolution of our observations (0.5”) im¬ 
plies a resolution of 625 AU at the distance to the Lagoon 
Nebula. Our observations are unable to rule out the presence 
of binaries having separations less than this distance. The ef¬ 
fect of undetected binaries in our sample may be understood 
if we consider the change in the properties of the star if it is 
physically an equal mass binary system. There is no colour 



log Age (yr) 



Figure 15. The distribution of masses (a), ages (b) and M acc (c) 
determined using different stellar evolutionary models. Solid blue 
lines are results estimated using the Siess et al. (2000) models, 
dashed red lines using the Bressan et al. (2012) models and dotted 
green lines using the Tognelli et al. (2011) models. 


change, while all its magnitudes are 0.75 mag brighter. This 
effect is similar to shifting the model tracks to account for 
differences between different stellar evolutionary models. 


However, estimating the number of possible binary con¬ 
taminants in our sample is difficult as previous studies of 
TTauri binaries have been inconclusive. Bouwman et al.l 
I20061 suggest that close binaries may lose their discs signif¬ 
icantly faster based on mid-infrared observations, meaning 
that most disc bearing CTTS are not generally found in close 
binaries. Daemgen et al. (2013) resolved the components of 
19 T Tauri binaries in Chameleon I. The authors found that 
M acc estimates only for the closest binaries (< 25 AU) were 
significantly affected. 


4.7 Completeness of our CTTS sample 

The completeness limit of the photometry is estimated from 
the peak magnitude of the luminosity function. We find that 
peak magnitudes are approximately 19.8, 18.8 and 19.4 mag 
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in r, i, and Ha respectively. Using the reddening and dis¬ 
tance corrections of Section 4.1, this translates to a limit of 
about 0.3 Mq in the r versus (r — i ) CMD. This limit can 
be seen in the histogram of derived masses (Fig. |10^ ). 

In reality, the completeness of the final sample of CTTS 
is also affected by the EWho selection criterion. The EWhq 
selection criterion is designed to exclude interloping chro- 
mospherically active late-type stars. As a result, the EWhq 
selection threshold is larger than the lower limit of EWhq 
in CTTS (see Sec. 3 and White & Basri 2003). To estimate 
the number of CTTS missed due to our adopted selection 
threshold, we must first estimate the total number of La¬ 
goon Nebula members in the mass range we observe, and 
then estimate how many of them are accreting, i.e. CTTS. 


The total number of stars in the mass range we observe 
(0.2-2.2 Mq) can be estimated from the initial mass func¬ 
tion, provided the slope, maximum stellar mass and total 
cluster mass of region is known. We adopt the results of Pris- 
inzano et al. (2005), where the authors find a power law in¬ 
dex of 1.22±0.17, and a total cluster mass of ~930 Mq. The 


maximum known stellar mass in the region is 20 Mq (Wei- 
dner & Kroupa 20061. Based on these results, we estimate 


that ~ 560, 260, and 40 Lagoon nebula members are present 
in the 0.2-0.4 Mq, 0.4-1 Mq, 1-2 Mq mass ranges respec¬ 
tively. Prisinzano et al. (2007) find that around 40 per cent 
of member stars in the Lagoon Nebula are actively accret¬ 
ing at a given time, i.e. are expected to be identifiable as 
CTTS. Combining these results with the estimated number 
of members, we can estimate the total number of CTTS in 
the Lagoon Nebula. Comparing this estimated number with 
our sample, we find that our sample is complete up to 60, 90 
and 50 per cent in the 0.2-0.4Mq, 0.4-1 Mq, 1-2 Mq ranges 
respectively. 


4.8 Contaminants in our CTTS sample 

The positions of the Ha identified sample of CTTS in the 
near-infrared and (u — g) vs. (g — r ) colour diagrams con¬ 
cur with the expected positions of CTTS. The M acc mea¬ 
sured from Ha agrees well with those measured using ri¬ 
band photometry. No stars in our sample have positions be¬ 
yond the lOMyr isochrones, with an overwhelming major¬ 
ity falling around the 1 Myr isochrone suggesting that back¬ 
ground/foreground contamination in our sample is small. 
This fact, along with the multiple signatures of disc accretion 
are highly suggestive that the vast majority of our candidate 
CTTS are genuinely accreting. Moreover, the EWh 0 crite¬ 
rion is designed to exclude the expected interlopers such as 
late type stars whose EWho is due to chromospheric activ¬ 
ity. However, we do not explicitly select for cluster member¬ 
ship (using for example proper motions). There are only four 
stars (out of 91) that are accreting in the (r — i) vs. (r—Ha) 
diagram and not accreting in the (u — g) vs. (g — r) diagram. 
Neglecting that such a change is due to accretion variability 
it is possible that these stars may be foreground unreddened 
late-type stars. Due to the nebulosity, background stars are 
not expected to be found (see Section 1). 

If we assume that the fraction of Ha emission stars not 
showing any u-band excess are field stars, then our contami¬ 
nation rate is ~ 5 per cent. The fact that these stars may not 
be accreting is strengthened by the fact that most of these 



Figure 16. The distribution of M acc with stellar mass (M„) of 
CTTS is shown by blue crosses. The best-fit power-law is given by 
Macc oc M* 2 ' 14 ^ 0 ' 3 . Red arrows give the lower detection limits. 


stars do not show atypical near-infrared excess. Note that 
there is no preferential spatial location for these stars. On 
the other hand, if we take the number of stars falling at the 
MS locus in the near-infrared to be contaminants, then our 
contamination rate ~15 per cent. However, a few stars that 
posses discs and are still accreting may not show infrared 
excesses due to variations in inclination angles (Meyer et 
al. 1997), or because they might be accreting from a debris 
disc. Moreover, it is difficult to differentiate CTTS using 
near-infrared colours alone, as the median (H — K ) colour 
disc excess from the UKIDSS survey in the near-infrared is 
only slightly larger than atypical main sequence colours (Lu¬ 
cas et al. 2008). In reality, our contamination rate is most 
likely somewhere in between ~ 5-15 per cent, and the to¬ 
tal number of contaminants is expected to be around 10-30 
per cent. 

Our estimations suggest that our sample is clean, and 
complete to a sufficient degree to estimate the statistical 
properties of the PMS stars within the Lagoon Nebula. 


5 Macc AS A FUNCTION OF STELLAR MASS 


The A/ acc rate is plotted as a function of M* in Fig. 16 We 


employed survival analysis linear regression using ASURV 


(Lavalley et al. 19921 to account for the lower detec¬ 


tion limits in fitting a power law function over the entire 
mass range. The lower limit of non-detections is calculated 
adopting the EWh» selection criterion for all points on 
an 1 Myr isochrone. We found the best-fit power law in¬ 
dex a = 2.14 ± 0.3. Not including the lower limits lead to 
slightly steeper slope. This strong dependence is well docu¬ 
mented in a few star forming regions over a wide range of 


masses (0.1-3 Mq; see Muzerolle et al. 2003 Natta et al. 
20061. The median a ~ 1.8-2.2 ( Hartmann et al.||2006 1 , al- 


though much steeper (|Fang et al. ||2009[ ) or shallower values 
( Barentsen et al.|201l| ) have also been found. 

It is important to consider the observational selection 


effects that may affect the observed relation (Clarke & 
|Pringle||2006| |Da Rio et al.|[2014| ) . The lack of T Tauri can¬ 
didates at the lower mass end in Fig. |16| is due to the data 
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reaching the detection limit, but we should be able to de¬ 
tect more CTTS at M* > 0.5 Mq with lower M acc - Stars 
with M* ~ 0.5 Mq are typically a magnitude brighter than 
the faint limit. Their paucity cannot be explained by the ob¬ 
servational selection limits. Similarly, it has been suggested 
that the upper envelope of accretion rates at any given mass 
is claimed to be due to the difficulty of detecting stars with 
L acc > L* (Clarke & Pringle 2006). However, such continuum 
stars constitute only a small fraction of the known popula¬ 
tion of PMS stars (Hartmann 2008). A small fraction of high 
M acc across the mass range will bias the observed relations, 
but cannot explain it. Therefore, observational limits can 
bias the observed relation, but cannot easily explain it. 


Overall, there exists a spread greater than an order of 
magnitude in accretion rates at any given mass in Fig. |16[ 
In Fig. |12| the accretion rates of 87 CTTS from the u and 
Ha excess emission have a mean difference <0.17 dex, much 
smaller than the spread in accretion rates. There are more 
likely other intrinsic factors besides stellar mass controlling 
accretion rates. Due to the ages of individual stars being 
uncertain, we do not examine the variation of M ac c as a 
function of individual stellar ages. 


SPATIAL DISTRIBUTION OF CTTS 
CANDIDATES 


Fig. |I7| shows the spatial distribution of CTTS candidates. 
Known O-type stars in the field are also shown. CTTS candi¬ 
dates are concentrated in multiple sub-clusters (having den¬ 
sity > 10 stars pc -2 ) in four regions: 

(i) in the southern bright rim M8 E. Tothill et al. (2002) 
identified several gas clumps within this region. A Herbig- 
Haro object has also been reportedly found by Arias et al. 
(2006). The 06 Ve star HD 165052, which is a confirmed 
member of the Lagoon nebula (Tothill et al. 2002) lies at 
the eastern edge of the rim. 

(ii) the embedded central cluster NGC 6530, which has a 
significant number of CTTS located to its immediate north. 


(iii) the Hourglass nebula region (marked as HG in 171, 
and a region to its north east. It contains at least two O-type 
members, the 04 V star 9Sgr, and the 07 V star H36 which 
are thought to be the chiefly responsible for ionising the La¬ 
goon and Hourglass nebula respectively. 9Sgr has 10 CTTS 
identified within a 1 pc radius. The number of CTTS sur¬ 
rounding H 36 is thought to be similar (see Arias et al. 2006), 
but cannot be confirmed in our study due to the difficulty 
in subtracting the highly variable background around H36. 
Both W hite et al.| (1997| and Tothill et al. (2002) report a 
significant amount of molecular gas (10-30 Mq) still shield¬ 
ing H 36, forming a compact HII region inside the molecular 
cloud. 


(iv) a north-western region (marked as NW in 171 at 
the edge of the visible nebulosity, surrounding the 09 star 
HD 164536, which is not a confirmed member. This region 
has received comparably little interest in previous studies of 
the Lagoon Nebula, owing to the brighter core, where star 
formation appears to be most active. 

A small fraction (< 6 %) of relatively more isolated stars 
are situated along the edges of the Lagoon Nebula, wherein 


some cases the nearest CTTS candidate is located at dis¬ 
tances > 1 pc. The overall morphology of our sample is in¬ 
dicative of a continuum ranging from dense clustering to 
comparative isolation. This is similar to results in nearby 
star-forming regions (Allen et al. 2007). Our spatial over 
densities of CTTS candidates are very similar to those found 
by Kuhn et al. (20141. 


6.1 Distribution of accretion rates 


It is interesting to see whether the spatial morphology of 
CTTS candidates is reflected in the distribution of accre¬ 
tion rates, and if any inference of how star formation has 
proceeded can be drawn from it. Lada et al. (1976) pro¬ 
posed star formation proceeded east to west through time 
beginning in NGC 6530 which triggered star formation in 
the Hourglass nebula. Lightfoot et al. (19841 put forth a 
similar scenario, where star formation was triggered by the 
NGC 6530 cluster in the Hourglass nebula region, and cite 
the presence of peripheral O stars in the Lagoon nebula as 
evidence for a younger population. Alternatively, Damiani 
et al. (2004); Prisinzano et al. (2005) and Prisinzano et al. 
(2012) suggested that star formation progressed from the 
north to the south, and the north-eastern CTTS were formed 
initially and have currently shredded their molecular gas, 
and star formation proceeded towards the west and south¬ 
west. Arias et al. (2007) reported that the ages of CTTS 
decreased from the north to the south. Tothill et al. (2008) 
consider these results, and suggest that star formation may 
have proceeded outwards from the core, evidenced by age 
gradients found by Damiani et al. (2004) and Arias et al. 
(2007), and the young sources around the dense molecular 
cloud cores at the edges of the nebula. 


To investigate such scenarios, a 2D histogram of ac¬ 
cretion rates covering the entire mass range, and restricted 
to the 0.35 Mq < M„ < 0.6 Mq are shown in Fig. |18| The 
colour is representative of the median accretion rate of all 
the stars situated in a 0.5 pc 2 box. The diffuse Ha contours, 
which outline the nebular emission are also shown. The dis¬ 
tribution of the accretion rates is considered to be due to 
the star formation history of the Lagoon nebula. Younger 
stars with higher accretion rates are thought to be located 
closer to their birth locations, while older accreting CTTS 
may have had time to dynamically evolve and move away 
from their natal molecular cloud which is thought to be re¬ 
flected by larger spacing between them, and lower accretion 
rates. However, the accretion rate also depends steeply on 
stellar mass (see Fig. |16| ), which is why we corroborate our 
map covering the entire observed stellar mass range with one 
restricted to 0.35- 0.6 Mq CTTS, which includes 135 candi¬ 
dates. Based on the best-fit from Fig. |16[ stellar mass affects 
the Mace by ~ 0.4 dex in this range, so any larger variations 
of Mace are likely reflective of dependences on stellar ages. 
Therefore, the distribution of accretion rates can be an in¬ 
dicator of the star formation history within the region, and 
also reflective of any intrinsic age spread. 

In Fig. |18[ we find that the median accretion rates are 
highest in the M8E and Hourglass region (—7.2Moyr“ ), 
and decrease non-uniformly towards the edges of the neb¬ 
ula. The accretion rate of the central embedded cluster 
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Figure 17. CTTS candidates overplotted on the Ho, shown as inverted grey scale. The plotted objects are coloured according to their 
estimated mass accretion rates. The positions of known O-type members are picked out by asteriks. The approximate centres of the 
Hourglass nebula (HG), NGC6530 cluster, the positions of the M8E rim and the north-western region (NW) are also labelled. The black 
bar marks an extent of 3.6 pc at the adopted distance to the Lagoon Nebula. 


NGC6530 is between —7.5 to —8.2MQyr _1 , and the north¬ 
western region fractionally lower ranging between —7.8 to 
—8.2 Mgyr -1 . 

Our results do not corroborate completely with any sin¬ 
gle sequential star formation scenario proposed in the liter¬ 
ature. If, we consider star formation to have proceeded east 
to west as described by Lada et al. (1976), why are the me¬ 
dian accretion rates in the south-eastern rim so high? Also, 
the region is known to contain at least one Herbig-Haro ob¬ 
ject (see Tothill et al. 2008), and dense CO clumps suggest¬ 
ing youth, if not younger then at least similar to the age 
of the Hourglass nebula region. Alternatively, if accretion 
proceeded north to south as suggested by Prisinzano et al. 
(2005; 2012), the median accretion rate of the north-western 
region is comparable to the central NGC 6530 cluster located 
to its south. The north-western region has not been as com¬ 
prehensively studied as the central Nebula, but |Kumar ~&1 
Anandarao] ( |2010| l identify six Class I objects using Spitzer 
infrared observations near the O-type star HD 164536, and 
comparable to the number density of Class I objects found 
in the NGC 6530 cluster. This gives weight to Fig. |18| that 
accretion in the north-eastern edge of the Lagoon Nebula 
is proceeding at rates fractionally lower to the central clus¬ 
ter. Alternatively, if star formation proceeded outwards from 
the core of NGC 6530 as suggested by Tothill et al. (2008) 


why are older stars found towards the northern, and east¬ 
ern edges, but not in the southern and western edges? It is 
possible that low M acc CTTS in regions of comparatively 
higher values of extinction are not detected in the study, 
which could lead to observed spatial dependencies. But the 
overall spread in measured extinction values within the re¬ 
gion is small (Ay ~ 0.5 mag, see Section 1). In our mass re¬ 
stricted M acc map, we find similar dependencies even though 
the imposed limit of 0.35 M@ places any star with extinction 
Ay < 2 mag within our detection limits. Both these consider¬ 
ations suggest that differential extinction is unlikely to bias 
our results significantly. Our results lead us to suggest that 
a simple sequential star formation process may be too sim¬ 
plistic to apply throughout the Lagoon Nebula, although 
there is evidence that some regions of the Lagoon nebula 
may be older than others. Such a conclusion was discussed 
briefly in Getman et al. (2014) in their multi-region study 
conducted using archive infrared and X-ray data. Instead, 
the distribution of accretion rates may be reflective of the 
clumped nature of the collapsing natal molecular cloud and 
the proximity of ionising stars. The reasons for cloud collapse 
(or triggering) of star formation may be to due to multiple 
factors, and therefore not explained by a single scenario. 
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Figure 18. Density distribution of M acc in the complete mass range (top) and in the restricted range 0.35 Mq < M* < 0.6 Mg (bottom), 
which contains 135 stars. The colour is representative of the median M acc within a 0.5 pc 2 box. White is regions with no identified CTTS. 
The corresponding colourbar is shown at the bottom. Known O-type stars (asteriks) are shown for reference. The approximate centres 
of the Hourglass nebula, NGC6530 cluster, and M8E region are annotated following Fig. |17| Overplotted are Ha image contours at 1%, 
5%, 10% of the total emission, which outline the diffuse nebular emission. 


6.2 Mapping of CTTS and protostellar 
populations 

The observed spatial distributions of young stellar objects 
(YSO) in star forming regions are a relic of the star forma¬ 
tion process, and also indicative of their dynamical evolu¬ 
tion from the protostellar to the pre-main sequence phase. 
Recent surveys of nearby star-forming regions (for example, 


see Winston et al.|2007 Megeath, Li & Nordlund 2009) have 


found young Class I protostars densely packed and tracing 
their natal molecular clouds, while relatively older pre-main 
sequence stars are found to be comparatively more spread 
out. Such morphologies place useful constraints on current 
star formation theories (Megeath et al. 2009) |Moeckel fe| 
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Figure 19. Distribution of CTTS candidates (circles), and Class I objects (magenta crosses) identified by Kumar Sz Anandrao (2010) 
overlaid on Spitzer 8pm contours. Colour of the circles is representative of measured log A fcci as in Fig.llTj The approximate centres of 
the four regions are annotated following Fig. EH 



Figure 20. Inverted Ho; greyscale of the M8E region overlaid by SCUBA 850pm contours. Symbols are same as Fig. |l9[ The O-type 
stars 9 Sgr and H 36 are located to the northwest. 


|Bonnell||2009| |Higuchi et al.|[2013| ), and can also be useful 
diagnostics of the evolutionary stage of PMS stars. 

The sample of Class I objects is drawn from a Spitzer 
study of the region by Kumar & Anandarao (2010). Briefly, 
the authors of that paper used observations in the Spitzer 
3.6, 4.5, 5.8 and 8 pm bands to classify PMS objects accord¬ 
ing to their spectral slope. They identified 65 Class I objects 
on this basis. The distribution of our VPHAS+ CTTS sam¬ 
ple, and Class I sources from Kumar & Anandarao (2010) 


are over-plotted on contours of Spitzer 8 pm emission in 
Fig. [19] The region surrounding the Hourglass nebula is 
saturated in the Spitzer image. Contours display a central 
egg shaped morphology, with elongated filamentary struc¬ 
tures found across the northern, southern and western edges. 
Class I objects are found within the filamentary gas struc¬ 
tures indicated by the contours. CTTS are found along the 
edges of the contour lines, and appear to hug the gas struc- 
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Figure 21. Inverted Ha greyscale image of an elephant’s trunk 
like structure in the Northwest of the Lagoon nebula region over¬ 
laid by Spitzer 8pm contours. Symbols are same as Fig. |19| The 
O-type star HD 164536 is located to the south west, and is possi¬ 
bly ionising the western rims of the structure. 


tures, with extended distributions in halos surrounding the 
Class I objects. 

The distribution of the nearest neighbours of Class I ob¬ 
jects has a median projected distance of 0.18 pc. The CTTS 
median nearest neighbour separation is 0.24 pc. There are 
no defined peaks, but an exponential decrease after a grad¬ 
ual peak. On performing a Kolmogorov-Smirnov (K-S) test, 
we find the probability that Class I and Class II objects are 
derived from the same random distribution to be 10 %, i.e. 
non-negligible (p-value = 0.09). We also generate a random 
distribution of 10,000 stars, where no underlying distribu¬ 
tion is assumed and we calculate their nearest neighbour 
separations. Comparing the YSO distribution with the sam¬ 
ple of random stars, the p-value in both cases is <0.01 %. 
The K-S tests suggest that the Class I objects and CTTS 
are similar. Both cannot be drawn from a randomly dis¬ 
tributed population of stars, suggesting significant cluster¬ 
ing. We caution that our numbers may be biased as both 
the CTTS and Class I samples are selected from data which 
are affected in the Hourglass nebula, thus possibly not iden¬ 
tifying important sources within this region. 


Adopting the mean molecular cloud properties of the 
south eastern gas clumps derived by Tothill et al. (2002), 
where the mean molecular density is 4.46 x 10 4 cm -3 , and 
mean gas temperature around 25 K we find the Jeans length 
( Jeans|1928[ ) to be 0.15 pc. The median observed spacing of 
Class I objects is roughly similar to the Jeans length. Similar 


results were found in nearby star-forming regions by (Young 
|et al.|[2006| |Winston et al.|[2007| |Enoch et al.||2009| ). How 
ever, a number of Class I objects are found within the median 
spacing. This has lead Winston et al. (2007) to suggest the 
possibility of Bondi-Hoyle or competitive accretion in pro¬ 
tostars in such regions due to the possibility of overlapping 
envelopes. 


Another important morphology to consider here is the 
location of the YSO with respect to the ionising stars. In 
Fig.|20[ 850 pm SCUBA contours spanning the southern fil¬ 
amentary arc (Tothill et al. 2002) are overlaid on the Ha im¬ 
age. Protostars are located within the molecular gas clumps, 


while CTTS are situated at the edges. The region is being 
ionised by the 04 V star 9 Sgr, and possibly the H 36 07 V 
star located inside the compact Hu region to the north west 
of the region. The position of the Class I objects and CTTS 
is suggestive of a decreasing age gradient north to south with 
respect to the central cluster NGC6530 and the Hourglass 
nebula. 


We identify an elephant’s trunk like structure in the 
north west (Fig. 211. CTTS are located at the outer edges 
of the structure. The structure may be being ionised by 
the O type star HD 164536 situated around 1 pc to its south 
west, causing the bright rims on the western side of the struc¬ 
ture. Most of the molecular gas in the structure has been 
evaporated, as suggested by the low intensity of Spitzer con¬ 
tours. Compared to the M8E feature, the CTTS here are 
located on the other side of the bright rim being ionised, 
which is a curious morphology. No Class I sources have been 
identified within the region by Kumar & Anandarao (2010). 
CTTS all have an age between 1.6 and 2.3 Myr. The median 
radial velocity of members within the region (0.6 pc Myr -1 ; 
Prisinzano et al. 2007), suggests that stars are located close 
to where they formed. The difference in the morphology of 
Class I YSO and CTTS between these two regions suggests 
that there is a lack of evidence to apply a global star forma¬ 
tion scenario throughout the Lagoon nebula. 


6.3 Accretion properties with respect to O stars 


PMS accretion is subject to the influence on its environ¬ 
ment by nearby ionising stars. Ionising radiation from O- 
type stars is thought to disperse nearby molecular gas, and 
remove circumstellar material from PMS accretion discs dis¬ 
rupting star formation. There has been evidence for the pho¬ 
toevaporation process decreasing the disc fraction of PMS 
stars within a few parsecs ( |Guarcello et al. 2010; Roccatagli-| 
ata et al.|2011| ). Alternatively, ionising fronts of O-stars are 
thought to trigger star formation by collapsing molecular 


clouds (see Elmegreen]2000 

Sicilia-Aguilar et al.|2010 

Bar- 

entsen et al. 2011, Karr, Manoj & Ohashi 2009; Snider et al. 

2009 

)■ 


To study the effects of O-type stars on the PMS popu¬ 
lation in the Lagoon nebula, we analyse the radial distribu¬ 
tions of the CTTS number density, accretion and disc prop¬ 
erties with respect to their distance from the known O-type 
members in Figs. [22] and |23| Two O-type stars are consid¬ 
ered, namely the 04 V type star 9 Sgr, which is thought to 
be ionizing the central part of the Lagoon Nebula, and re¬ 
sponsible for clearing gas in the rift within the NGC 6530 
cluster (Tothill et al. 2008), and the 07 V type star H36 
which is thought to be the main ionizing force in the Hour¬ 
glass nebula (Arias et al. 2007). 


Fig. [22] indicates that the number density of CTTS 
within 1 pc of 9 Sgr is considerably higher, with the Spear¬ 
mans p=— 0.72. The median EWhh and logM acc are also 
considerably different from the mean within 1 pc of 9 Sgr 
(where pis 0.68 and -0.84 respectively), and follow a steady 
distribution 2cr within the median thereafter. We find little 
change in the median (H — K) value with radial distance 
from 9 Sgr except at distances >5 pc, at the edge of the neb¬ 
ula. The CTTS disc fraction (defined as the ratio of CTTS 


falling within the near-infrared excess region of (Lada et al. 
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Figure 22. Histogram showing the radial variation of number, number density, median EWh 0 , median log Af acc , median {H — K) colour, 
and disc fraction with respect to the O-type star 9Sgr (a-f). 
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Figure 23. Histogram showing the radial variation of number, number density, median EWh q , median log Af acc , median {II — K) colour, 
and disc fraction with respect to the O-type star H36 (a-f). 


) to the whole CTTS population within the region) fol- 
a similar trend, with the disc fraction somewhat smaller 
within 1 pc but steady thereafter until beyond 5 pc. 

We find that the number density is fractionally higher 
within 1 pc of H 36 (Fig. |23| ) , although this region is affected 
by incorrect Ha background sky subtractions (see Section 
2). If the number of known CTTS (from Arias et al. 2007) 
within the excluded region is considered, the number density 
increases to 72 %. We find the median median EWho and 
log M acc are higher within 1 pc. The disc fraction within 1 pc 
is significantly smaller, although we caution that our obser¬ 
vations are biased in this region (see Section 2.1). Arias et 
al. (2006) used near-infrared observations to find that the 


fraction of disc-bearing members to non-members is around 
70 % within this region, suggesting that it is close to the 
median in the Lagoon Nebula. 

Overall, we find that the gradients of CTTS accretion 
properties with respect to the O-type stars are not indica¬ 
tive of triggered star formation. We find evidence for disc 
photoevaporation within 1 pc of 9 Sgr and H 36, although 
we caution that taking into account censored observations, 
this may not be the case around H36. Finally, we find that 
the number density, median EWho and logM acc are com¬ 
paratively higher within 1 pc of both 9 Sgr and H36. This 
may be explained by the fact that the youngest stars are 
situated close to the O-type stars, as shown by Fig 118| 















22 V. M. Kalari et al. 


7 SUMMARY 

We identify 235 Classical T Tauri candidate stars in the La¬ 
goon Nebula and measured their accretion rates. The ma¬ 
jority of our sample demonstrate multiple indicators of disc 
accretion: Ha emission, UV-ex cess, and near-infrared ex¬ 
cesses. The conclusions drawn from our study are: 

(i) A majority of the CTTS in the Lagoon Nebula have 
isochronal ages around 1 Myr. This is consistent with the 
main-sequence lifetimes of O-type star members, the dynam¬ 
ical evolutionary stage, and the lifetime of the HII region. 

(ii) The accretion rates determined from Ha line luminos¬ 
ity and UV-excess luminosity from 87 stars show no statis¬ 
tically significant differences. The mean amplitude variation 
~ 0.17 dex. 

(iii) We find best-fit power law relation M acc oc A/„ 2 ' 1±0 ' 3 . 

(iv) The distribution of CTTS is a continuum ranging 
from dense clustering to relative isolation. Class I objects 
are comparatively more clustered than CTTS. 

(v) The youngest accretors are concentrated around the 
Hourglass nebula and M8 E region. Accretion rates decrease 
towards the edges of the nebula, except in the northwest. 

(vi) We find no evidence supporting a global star forma¬ 
tion scenario within the Lagoon nebula. We suggest that 
star formation may be influenced by environmental effects 
on smaller scales. 

Overall, we demonstrate the application of VPHAS+ data 
to identify and measure accretion rates of CTTS candidates 
using uniform selection criteria and known detection limits. 
Future astrometric and spectroscopic observations such as 
those undertaken by the GAIA and GES survey will obtain 
further information about a significant number of our CTTS 
sample, complementing our VPHAS+ study. 
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